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Graphity models are characterized by configuration spaces in which states corre- 
spond to graphs and Hamiltonians that depend on local properties of graphs such as 
the degrees of vertices and numbers of short cycles. As statistical systems, graphity 
models can be studied analytically by estimating their partition functions or numer- 
ically by Monte Carlo simulations. Results presented here are based on both of these 
approaches and give new information about the high- and low-temperature behavior 
of the models and the transitions between them. In particular, it is shown that mat- 
ter degrees of freedom must play an important role in order for the low-temperature 
regime to be described by graphs resembling interesting extended geometries. 



INTRODUCTION 



The idea of describing space (or spacetime) as a product of a dynamical process involving 
some discretized elements resonates with many approaches to stu dying th e unification of quan- 
tum mechanics with gravity [H, i, i, i, H, 0, 0,1, B 0, [HI, E, 0, H, H [13 • "GraphitV' 



refers to a class of models in which the degrees of freedom represent dynamical graphs 

mm 

and whose Hamiltonians depend only on minimal information encoded in graphs such as the 
degree of the vertices and cycles formed by the edges. It has been argued that such models 
can be useful in studying how an extended and/or latticelike geometry might emerge from a 
background independent dynamical system. 

Several approaches to studying quantum geometries are based on general or directed graphs 
and graphity models share some features with them while differing in others. In contrast 
to some setups defined primarily by a set of dynamical rules 0, Hf, 0, 0], graphity models 
have Hamiltonians, which allow them to be interpreted as statistical mechanical systems. In 
contrast with approaches starting from a quantization ansatz for general relativity jl, 0, [l^. 



111 . |12| . |13| , graphity models are only based on graph-theoretic notions and are not specific to 
a particular dimension or embedding manifold. They are in nature most similar to dynamical 
triangulations js], [o], 13] in the sense that their macroscopic properties are determined by the 



collective interactions between a large number of small elements. The graph models are also 
amenable to simulation just like dynamical triangulations js], lo], However, since not all 
general graphs correspond to triangulations, graphity models have a larger configuration space 
and their defining Hamiltonian cannot be easily mapped to a discretized version of the action 
of general relativity. 

The motivation for studying such general systems comes from asking to what extent geomet- 
ric notions such as the continuum, dimension, or macroscopic locality, should be considered as 
fundamental (as they are in all approaches that postulate the existence of a space or space-time 
manifold) and to what extent they may be only approximate or emergent. One way of probing 
these issues is to consider a setup in which none of these notions is obviously true. Graphity 
models are an attempt in this direction that assume only the existence of a Hamiltonian function 
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H determined by micro-local conditions of graphs (for a discussion of the distinction between 
macro and micro locahty, see jl^.) 

It has been argued that graphity models can behave very differently at high and low energies 



17|, [ISj. At low energies, the preferred graph state can represent an emergent geometry, while 
at high energies, the system is expected to be disordered and highly connected. The purpose 
of this paper is to study the statistical mechanics of the models in order to learn more about 
the nature of the low-temperature regime. The main finding is that while the most simple 
models do generate extended graphs at low temperatures, these graphs are effectively only one 
dimensional. However, when matter degrees of freedom that are sensitive to graph boundary 
conditions are added, the low-energy behavior can dramatically change and various homogenous 
latticeli ke g raphs r nay appear. It thus appears that matter, as has been noticed in other contexts 
as well [la 20, 21, 22], should be regarded as an important, if not crucial, factor in discussions 
of emergent geometries. 

The definition of graphity models is reviewed in the next section. A distinction is made 
between a basic, or bare-bones, model of pure graphs and more sophisticated versions that can 
incorporate matter degrees of freedom in the form of gauge fields as well. The interpretation of 
the models as statistical systems is also introduced laying the ground for much of the following 
discussion. Section IIIII deals with the partition function of the basic model. This discussion 
is detailed and involves identifying the types of graphs that can contribute significantly to the 
system's partition function. The discussion should be seen as complementary to the results of 
Monte Carlo simulations presented in Sec. IIV[ Indeed, it is the simulations that motivate and 
at the same time justify the study of particular types and the omission of other types of graphs 
in the theoretical section. As such, the analysis in Sec. IIIII is an attempt to provide analytical 
understanding of the results of the numerical simulations. And in turn the analytical work 
suggests what kind of numerical studies can be of potential interest. 

The results of the Monte Carlo simulations in Sec. IIVI indicate that the basic graphity 
model is insufficient to generate extended geometries with dimension greater than one in the 
low-energy regime. The emerging graphs are either chainlike or treelike. It is shown in this 
section, however, that supplementing the basic model with a condition for graph homogeneity 
does lead to very different and interesting results. This observation motivates the study of graph 
models with matter in Sec. |Vl The discussion there mirrors that in Sec. IIIII and demonstrates 
that under certain quite general assumptions regarding the energy spectrum of the matter 
content, a model with a realistic matter component may indeed produce extended geometries 
at low temperature. A summary of all the findings is presented in Sec. I VI I 



II. GRAPHITY MODELS 



Graphity [17|, |18| is a name for a class of models in which states are associated with graphs 
and the configuration space is the space of all possible simple graphs with a given number of 
nodes A^. The models are defined in a quantum mechanical context so that the configuration 
space is actually a Hilbert space. Indeed the Hamiltonian of the model is most naturally 
formulated as a quantum mechanical operator rather than as a classical object. For this reason, 
the models can also be referred to as "quantum" graphity. 
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The total Hilbert space for a graphity model can be decomposed as 

N{N-l)/2 N 

'Htotal = 'Hedge ^ vertex ■ \^ ) 

The tensor products of TC^ertex and Hedge corresponds to putting degrees of freedom on each 
vertex and every possible edge connecting the vertices. Specific graphity models are defined by 
particular choices for the Ti-vertex and Tiedgei along with an appropriate Hamiltonian acting on 
the resulting Htotai- 



Basic Model 



In the basic model described in [18| , the Hilbert space associated with the vertices is trivial 



and the Hilbert space associated with each edge is spanned by two vectors 

7^ed3e = span{|0), |1)}. (2) 
The total Hilbert space is therefore 

N{N~l)/2 

Htotai = Hedge- (3) 

Denoting states in each copy of Hedge as \nab) with Uab = 0, 1, the total Hilbert space can be 
said to be spanned by the product states 

Htotai = span{|ni2) ® Iriis) (g) 1^23) ® ■ ■ ■ } (4) 

A general state in this space is a superposition of the basis states with complex coefficients. 

Interpreting the states |1) and |0) as conveying whether or not, respectively, a link is 
present between the two vertices, each of the basis states (jl]) can be associated with a graph 
configuration or a graph diagram. Thus, the total Hilbert space (13]) can be decomposed as 

Hedges = ^ Hq (5) 
G 

where He is a space spanned by a single basis state and corresponds to a single graph G. 
Because the vertices are labeled and are thus distinguishable, the tensor summation over Hg 
can contain multiple graphs G that are isomorphic to each other. 

As argued in |18j, a Hermitian Hamiltonian H acting on Htotai can be used to associate an 
energy E{G) with a graph states \iPg) through the relation 

E{G) = {^g\:H: I^g). (6) 

The form of the Hamiltonian is restricted by requirements of general graph locality and can 



thus depend only on a few properties of the graph. The basic model described in [18 
Hamiltonian of the form 



uses a 



H = Hy + Hb + Hint (7) 
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where Hy depends on the degree of vertices, Hb depends on the cycle structure of a graph, 
and Hint is an interaction term. Each of the terms in the Hamiltonian can be implemented 
using creation and annihilation operators acting on Htotai [S]- For the purposes of this paper, 
however, it is sufficient to describe only the effective properties of the energies E associated 
with these terms. 

The valence term Hy assigns an energy to the graph according to 



Ev =9vY^ 



Here gv is a positive coupling constant, p is some positive real number, v{a) is the degree of 
vertex a, and Vq is positive integer that determines the preferred valence of each vertex in the 
graph. 

The other term Hb is such that its contribution to the energy can be written as 

EB = J2EBia) (9) 

a 

with 

N{N-l)/2 

EB{a) = - 9B{L)P{a,L). (10) 

L=3 

The last expression is a contribution that corresponds to a single vertex a and that depends on 
the number of cycles P{a, L) of length L that pass through that vertex. A cycle is defined as a 
closed walk along the edges of a graph in which no edge is traversed more than once (see [3] 
for more details.) The lower and upper limits on the sum are a consequence of this particular 
definition of a cycle. 

The effective coupling gB{L) associated with cycles of each length is 

9B{L)=gB'^ (11) 

with gB and r some constants. The form of this effective coupling plays a very important role 
in the properties of graphity models. The interplay between the exponential in the numerator 
and the factorial in the denominator of (fTT!) means that the effective coupling increases with 
L for small L but decreases very rapidly for large L. The rapid falloff implies that the energy 
Eb is primarily determined by the number of short cycles in a graph. The total energy (ITU]) 
may be well approximated by restricting the summation over L to the range <C L <c Lmax 
follows: 

Lmax 

EB{a)^-Y,gB{L)P{a,L). (12) 

The maximal cycle length Lmax niust be determined for each model depending on the given 
parameters vq and r and on what accuracy the energy is to be evaluated. 

The turnover between increasing and decreasing gB{L) occurs at some length L^, determined 
by the parameter r. When the effective coupling gB{L) is multiplied by the number of cycles 
P{a,L), which usually increases with L, the maximum value of the product occurs at another 



(c) 

FIG. 1: Interaction moves on graphs. 



length L^^ determined by both r and the type of graph used to extract P(a, L). Both scales 
and L^K^K are smaller than Lmax- It is important that none of these cycle length scales depend 
on the number of vertices in the graph A^. Therefore, since only short cycles are important, the 
energy function can be said to be quasi- (micro) local. 

Other important properties of (fTTl) are determined by the signs of r and qb- For r > and 
qb > ^1 the effective coupling is positive definite; when it is inserted into the energy formula 
it is found that all cycles contribute a negative energy. When r < and > 0, the sign 
of gB^L) is different for even and odd cycles; even cycles contribute negative energy, and odd 
cycles contribute positive energy. Cases with qb < can be understood similarly. 

Moving on from the definition of the cycle term, note that the basis states (jl]) are all 
eigenstates of the operators Hy and Hb- In order for the system to evolve from one graph 
configuration to another, an interaction term is necessary. While many forms of interactions 
acting on graph states are possible, the ones suggested in [l7|, [l8| are among those that act 
locally on graphs. Examples of such interactions are shown in Fig. [TJ The first two have 
the property that they preserve the degree of all vertices; the third move allows the degree to 
change. All of them preserve the connectedness of a graph state. 



B. Model with Matter 



Extensions of the basic model may involve introducing a nontrivial Hilbert space Tivertices 
for the vertices or exten ding th e Hilbert space Tiedge associated with the edges. An example of 



the latter, described in [17|, ll8|, consists of splitting the |1) state in the basic model into three 
distinct states so that 

Tied^e = span{|0), | li), | I2), | I3)}. (13) 

The states denote "on" edges with different internal labels. 

A Hamiltonian for the extended model can depend on the internal degrees of freedom of 
the on edges. Thus, in addition to analogs of the terms described in ([7]), other terms are also 
possible giving a total Hamiltonian of the from 

H = Hv + Hb + H^nt + Hm + Hg. (14) 

Here the notation Hm and Hq is meant to suggest that these terms may give rise to effective 
matter {Hm) and perhaps gravitational {Hg) dynamics. These terms could be implemented. 
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as originally suggested in 17, 3], using the string-net mechanism 23, 24, 25 1 related to the 
Kogut-Susskind formulation of gauge theory j26|, or perhaps by other means, for example, by 
adapting the framework of algebraic quantum gravity to dynamical networks [l3| . 

The specifics of the implementations of these terms will not be important for the work 
in the present paper. However, it might be expected that a realistic implementation should 
contain wavelike excitations such as photons or quantum mechanical particle wavefunctions in 
the matter spectrum. This rather weak assumption will be invoked in Sec. |Vl 



C. Statistical Interpretation 

As a statistical system, a graphity model can be studied by computing its partition function 
in the canonical ensemble at a characteristic temperature T. As the configurations of the basic 
graphity model correspond to graphs, the partition function can be written as 



Z = Y, e-^^(^\ (15) 



G 

where 13 = l/ksT is the inverse temperature, fcs is the Boltzmann constant, and E{G) is the 
graph energy computed using ([6]). 

In models containing extra degrees of freedom M, perhaps corresponding to matter, the 
partition function must involve another summation of those states, 

Z = e-PE{G)-pEM{GM)_ (^X6) 

G M 

The energies Em of the additional degrees of freedom can depend on the underlying graph 
configuration G, and so this dependence is explicitly shown in f|T6|) . When the main subject 
of study are graphs rather than the matter degrees of freedom, the summation over M can be 
hidden by writing 

Z = ^ z(G)e-^^(«) (17) 

G 

with some function z{G). It will be useful to note that both f[T51) and flTTl) can be cast into the 
general form 

a 

if the quantity C{G) is defined as 

C{G) = PE{G)-hiz{G). (19) 



This quantity appearing in the exponent in the Boltzmann factor can be called the Boltzmann 
criterion. It reduces to the usual (3E{G) if the partition function for matter z{G) is set to unity. 
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III. STATISTICAL MECHANICS OF THE BASIC MODEL 

The purpose of this section is to estimate the partition function of a graphity model and 
use it to understand the statistical mechanics and "phase" structure of the model. Attention 
is restricted to classical graphity models so that the sum in the partition function is over the 
basis vectors of Titotai only, i.e. states that correspond to classical graphs. Furthermore, the 
matter factor is set to 

z{G) = 1 (20) 

for simplicity, effectively putting the partition function in the basic form (fT5|) . 

Since the aim of a graphity model is to describe an emergent geometry, emphasis is placed 
on understanding the partition function in the limit of a large number of vertices, N ^ oo. 
Since the models considered all have a finite number of vertices, this limit should be thought 
of in the following way. Consider a family of models each with a different but with the same 
other basic parameters (couplings gv, Qb or parameters vq, r). The equilibrium properties of 
each of these models can be evaluated separately. The continuum limit is understood as the 
limit of this series of systems arranged in order of increasing A^. 

The estimate for the partition function below is very crude and is based on only a few 
types of graphs. These are described in Sec. IIII A[ It is argued that these few configurations 
are sufficient to extract basic information of the model such as the approximate transition 
temperatures. These properties, and the limit of applicability of the results, are described in 
Sec. UTTBl 

A. Configurations 

In any statistical system, some configurations have more importance than others in the 
partition function and thus in determining the expectation values of observables. In the context 
of graph models this means that some graphs, or classes of graphs, contribute more than others 
to Z. Understanding a graphity model can therefore be phrased as the problem of identifying 
the dominant classes of graphs and evaluating their contributions to Z. 

In this section, it is argued that the important classes of graphs to consider are random 
graphs and graphs with lowest energy. Other types of graphs that are of interest and are also 
discussed are latticelike graphs. Thus, the partition function is split into 

Z = Zn + ZL + ZH + --- (21) 

Each of these terms represents a sum like (fTSll restricted to a special class of graphs and will 
be estimated separately. The ellipsis denotes contributions from other graphs that do not fit 
into the mentioned categories. 

In what follows the components of Z in 0211) are written for large but finite A^. This means 
that often only the dominant behavior of each term, as a function of and f3, is presented. 
Unless otherwise stated, estimates for InZ are valid up to terms or order A^, which can be 
neglected given that the dominant contributions are usually of order A^lnA^. Dependence on 
the inverse temperature (3 is always shown explicitly even when it appears subleading as a 
function of A^. 
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1. Low- Energy Graphs 

A class of graphs that can contribute significantly to the partition function are those which 
have very small energy and thus a very large Boltzmann factor. Since the energy of a graph is 
determined by the valence and the cycle structure, minimizing the energy must take both of 
these properties into account. Consider regular graphs with varying degrees v. The number of 
cycles P(a, L) of length L at each vertex a is bounded from above by 

P{a,L)<cv^-^ (22) 

for some constant c. Since the bound is generous, the value of this constant is not important 
and can be set to unity. It follows that the total energy due to cycles is bounded by 

E 9b{L) Pia, L)<gBN . (23) 

a L 

The absolute value in the exponent is required for the expression to hold for both negative and 
positive r. Putting this bound together with the energy contribution from vertex degrees, one 
obtains another bound £'o,^ for the total energy of regular graphs of degree v 

E{G) > Eo,, = ^y7Vef(--o)' - geN^. (24) 

The first term is minimum for v = Vq and grows superexponentially for v different than Vq. The 
magnitude of the second term increases with v but does so more slowly than the first term. 
Thus, it is possible to arrange the parameters so that the lowest possible energy occurs when 
the degree is vq. To do this, note that the energy of a I'o-regular graph is smaller than 

E,, = gvN (25) 

and the energy of a + 1-regular graph must be larger than 

ekl("o+i) 

E,,+i = gvNeP - geN — -. (26) 

Setting EyQ^i — E^^ > one finds 

e\r\{vo+l) 

gv {eP - 1) - gs — - > 0. (27) 

vo + l 

This is a conservative requirement on the couplings gv and gs that guarantees the lowest 
possible energy to be associated with a Vo-regular graph (as opposed to a v- regular graph). 
When the inequality is satisfied, it also implies that the smallest energy graph is exactly regular - 
all graphs which break the regularity condition even at a select number of vertices have energy 
higher than the minimal energy regular graph. The condition also ensures that the energy 
function is bounded from below in the N ^ oo limit. 
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In what follows, it will always be assumed that the couplings gs and gy are such that the 
minimal energy graph Gq is Wo-regular. The contribution of this graph to the partition function 
may be estimated as 

Zz.-e^'-^^e-^^^^"). (28) 
It will be convenient to write this and all the following estimates in the form 

IuZl N\nN - l3E{Go). (29) 

Here, E{Gq) is the minimum value of the energy associated to Gq. The factor involving A^lnA^ 
is due to the fact that there are on the order of N\ ways of assigning labels to the vertices 
of a graph with vertices; although these different assignments are sometimes considered as 
equivalent from the graph-theoretic perspective, they correspond to distinct states in Htot and 
thus must be counted in Zl- 

The energy of a graph E{Go) can also be written as 

E{Go) = N{gv + eo) (30) 

giving eo the interpretation as the average cycle energy per vertex in Gq. Because of the form 
of the cycle term in the Hamiltonian, this quantity can be (and in cases of interest, is) negative. 
With this notation, becomes 

\nZL^ N\nN -f]N{gv + eo). (31) 

There may be several other graphs G that have an energy E{G) very close to E{Gq). Such 
graphs could be interpreted as perturbations of Go and would contribute a quantity close to 
(l3T|l to the full partition function. This could be taken into account by multiplying Zl by the 
number of such perturbation graphs. However, this multiplicative factor should be expected to 
be on the order of A^^ for some value x that does not strongly depend on or /3. Thus, for 
large N, the effect of such a multiplicative factor will be negligible compared with the other 
factors in (1311) and hence will not be included in the analysis. 



2. Random Graphs 

Random graphs are configurations formed by randomly assigning edges to a set of vertices 
27(1 . Regular random graphs are those random graphs that are subject to the constraint that 
all the vertices have equal degree. They are widely studied in graph theory, and a large number 
of their properties are known. For example, a typical random graph with a large number of 
vertices and low degree has almost no short cycles j27[. Thus, a typical random graph G^^^ 
with degree v is assigned almost zero energy by the Hamiltonian Hb so that the total energy 
is entirely determined by the valence term 

EiGn,.) gvNe^^''-'"^^". (32) 

Thus, compared to graphs with low energy discussed previously, the Boltzmann factor of ran- 
dom graphs can be small. However, because there are many random graphs in the space of all 
possible regular graphs, they may contribute significantly to the partition function overall. 
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The number of regular random graphs is difficult to evaluate for general v and although 



some asymptotic formulae are known [27[. For the present purpose, an upper bound W{N,v) 
for the number of random graphs will be sufficient. In a f- regular graph with vertices, there 
are in total N{N — l)/2 edges of which Nv/2 are "on." The upper bound W{N,v) may be 
estimated by the binomial coefficient 

'^''('^.")=(iV(r-l)/2) (») 

This estimate overcounts graphs, because it does not take into account the requirement that the 
graph must be regular. When = 6 and f = 3, for example, this estimates gives PF(6, 3) ~ 10^"^ 
whereas there are only two nonisomorphic three-regular graphs with six vertices, for a total of 
2 ■ 6! = 1440 graphs. Nevertheless, the estimate becomes more accurate when is large. In 
that regime, W{N, v) can be expressed in a simpler form by writing out the binomial coefficient 
in terms of factorials and using the Stirling approximation. One finds 

\nW{N,v) ^]^NvhiN + 0{N). (34) 

The random graphs that minimize fl32l) are those with v = Vq, and those will be of primary 
interest here. Thus, using the observation fl5^ and the upper bound flMl) . the estimate for the 
contribution to Z from fo-regular random graphs is 

lnZ^,,„ ~ \nW{N,vo)-l3E{GR^J 

1 (35) 
= -NvolnN-pNgv. 

Within the random graph category, there are also graphs which contain more than Nvq/2 
edges. These graphs necessarily have at least two vertices with degree different than vq. Such 
graphs have minimum energy if they can be thought of as consisting of a f g-regular random 
graph with an extra set of edges. The number of such graphs could be estimated by taking 
W{N, vq) and multiplying by A^^^, where x is the number of additional edges. For small x, the 
energy of such a graph would be different from E{Gr^vo) xgvi^eP — 1). Thus, the contribution 
to the partition function of such graphs would be written as 

InZfi^^o,^ ~ ln^ij,^,o + 2xlnA^ - x(3gv{e^ - 1). (36) 

This formula can only be a reasonable estimate for small x. When x becomes comparable to 
A^, some vertices would have degree grater than fo + 1, and the estimated energy shift used in 
( !36|) would break down. In fact, the energy penalty for those graphs would grow exponentially, 
and hence they would contribute negligibly to Z. 

Summarizing the discussion on random graphs, their total contribution to Z is 



+ + (37) 



with X ranging from 1 to some number of order A^ and the ellipsis denoting further contributions 
from random graphs with more edges than Nvq/2 + x. These contributions can be a priori very 
significant, but since the goal of the graphity model is to describe graphs with low degree, the 
discussion below in Sec. IIIIBl is devoted to finding conditions under which these contributions 
can be neglected. The terms shown in (1371) will be sufficient for this purpose. 
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3. Homogeneous Graphs 



Lattice graphs are graphs that are regular and further are highly homogeneous. Because of 
this, these graphs are of great interest from the perspective of using a graphity model to describe 
an emergent geometry. However, it is important to note that latticelike graphs may or may 
not correspond to graphs with lowest energy. As such, the contribution Zh of homogeneous 
latticelike graphs should be considered separately from Zl- 

The estimated partition function Zh may be written as 



This has the same form as fl3Tl) but the additional sum is restricted (as denoted by the prime) 
to homogeneous, latticelike graphs G each of which has a cycle energy per vertex denoted by 
eniG). The number of terms in the restricted sum that are large depends on parameters Vq and 
r, but in general can be taken to be proportional to or some power thereof. In comparison 
with the dependence on N of the individual terms, however, this scaling with N is insignificant. 
Furthermore, since the exponents in the Boltzmann factor of (138|) are of the same form as in 
(l3T|l but the energy per node en is at most as small as eo, the homogeneous graph contribution 
will in general be much smaller than Zl for large N. 



The Hilbert space Htotai also supports a vast number of graphs that do not fit into the 
categories described above. Such graphs could be vq regular with some short cycles but not as 
many as Go, or almost regular with a few vertices breaking the valence condition, etc. They can 
be seen as configurations bridging the gaps between the random, low-energy, and homogeneous 
graphs. 

Writing down a partition function for these graphs is as difficult as writing one down for the 
full model and will not be attempted here. However, in certain regimes the contribution of these 
intermediate graphs can be ignored. In situations where random graphs dominate the partition 
function, the intermediate graphs are subdominant because they are fewer in number. And in 
situations where low-energy graphs dominate, the intermediate graphs are also subdominant, 
because they have higher energy. In Sec. IIVI numerical simulations of the full system will 
provide a posteriori justification for ignoring these intermediate graphs. 



Given the arguments in Sec. IIII Al the partition function of a graphity model may be 
approximated by the contribution due to random graphs and low-energy graphs. 





G 



4. Other Graphs 



B. Scaling and Thermodynamics 



Z^Zr + Zl. 



(39) 
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Using fl5^ . one can apply standard formulae 

U = (40) 

to find the average energy U or the specific heat C. Since Z has a number of terms that may 
dominate in some temperature range, one can expect a number of maxima in the specific heat 
that correspond to transitions between the dominance regions of these terms. The canonical 
way of finding the temperatures of these transitions would be to consider dC /dp = and solve 
for /3. However, the resulting equations are very involved, and the solutions can only be written 
neatly after using some approximations. So, instead, one can treat the problem in a different 
manner. 

Consider using only pairs of terms from fl39|) at a time to calculate the maxima in the specific 
heat. In particular, consider the pairs 

Zl,vo = Zl + Zr^v^, (42a) 
Zl,x = Zl + (42b) 
Zvo,x = Zr^vq + Zn^y^^i. (42c) 

Each of these partition functions has only two terms. Thus, specific heats derived from them 
have one maximum each marking transitions between dominance regimes of the two terms. For 
example, the function Zl^^^ has a characteristic inverse temperature that describes a transition 
from a fo-regular random graph regime (low j3) to a low-energy graph regime (large j3). As 
another example, the maximum in a specific heat derived from Z^^^^ describes a transition 
between graphs that should be expected to be Uo-regular and graphs that have x extra edges. 

In the limit N oo^ the maxima in the specific heats corresponding to fl42p diverge. The 
transition inverse temperatures are given, in the large limit, by 

/3l,.o = -i(^-l)lnAr (43a) 
eQ \ Z / 

Plx = ^ (l-l + '^hiN (43b) 

In the latter two equations, the quantity 

Aev = gv{e''-1) (44) 

is always positive and can be made as large as desired by adjusting gy- The quantity eo, being 
the energy per vertex coming from the cycle term, is proportional to the coupling qb, and in 
cases of interest is less than zero. At this stage, one can also note that the same transition 
temperatures could also be obtained by simply comparing the expressions for the components 
Zl and Z^i without computing the specific heats and their maxima. 
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As all the expressions contain a factor In A^, the transitions they describe all technically 
occur at zero temperature in the N —>■ oo limit. This is perhaps not surprising as it can be 
already seen in expressions like ([35D and flHT]) where the leading terms proportional to (3 are 
linear in N, while the leading terms independent of (3 scale as A^lniV. In the following, it is 
useful to either compare the solutions P at finite A^, or to divide out by InA^ to consider only 
the coefficients. 

In order to reconstruct the transitions in the full system, one can compare the solutions (143|) 
pairwise. Consider first the pair and Pvo,x- If the coupling condition (1271) holds (as is 
necessary for the estimate to be consistent), then 

I3l,vo > Pvo,x- (45) 

This implies that the lowest temperature regime is characterized by low-energy graphs and that 
there is a transition to exactly Wo-regular random graphs at an inverse temperature f3L,vo- The 
next transition at a much lower inverse temperature occurs between fQ-regular random 
graphs and random graphs that are not-quite regular. An important thing to note is that these 
transitions can be well separated when gy is large. Also, note that Py^^^ does not depend on x. 

Next, consider the pair Pl^^^ and Pl,x- When x is order unity (as is assumed) and A^ is 
large, these two quantities become equal. It would be incorrect, however, to conclude that 
the transition from low-energy graphs to not-quite regular random graphs occurs at the same 
temperature as the transition to exact random graphs. This is because solution P^^x already 
implies that the transition to not-quite regular graphs occurs at much lower values of p. The 
apparent paradox is an artifact of using the simplified function Z^^^, which does not take 
into account regular random graphs that form an intermediary stage between irregular graphs 
and low-energy graphs. The robust conclusion to draw from comparing Pl,vo and Pl,x is that 
low-energy graphs are stable at P higher than both of these values. 

Comparing the final pair of solutions and observing that Pl^x > Pvo,x would suggest that 
graphs are fo-regular at the highest temperature, irregular at intermediary temperature, and 
low-energy at low temperature. While this is consistent with the hierarchy of the solutions 
(H3|) . it is in confiict with the previous interpretation. Deciding between the two interpretations 
requires some extra input, and this can be obtained by calculating the expectation values of 
the energy around the transition Py^^x- It turns out that the original interpretation of the 
transitions with the irregular graphs as the highest energy stage is the correct one and, again, 
the apparent contradiction can be traced to the inappropriate use of the expression for Z^ x- 

The transition inverse temperatures fj^3|) must be understood as only approximations to 
those one would find by considering the better approximation fl39|) or the full graphity partition 
function. Nonetheless, the property that the transition /5s scale logarithmically with A^ should 
be consistent with what one would find from the full partition function. The estimate Py^^i for 
the lowest temperature transition should also be accurate for large A^ (as will be confirmed in 
the next section). And furthermore, the separation of the transitions from low-energy graphs to 
regular random graphs and from regular random graphs to irregular graphs should also persist 
in the full model. 

More specific questions regarding the nature of the transitions, however, cannot be addressed 
reliably. In particular, while (143|) are all associated with diverging specific heats, the presented 
analysis does not make it clear whether these divergences persist in the full model and whether 
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they can actually be associated with phase transitions. One possibility is that the lowest-energy 
transition described by Pl^^^ is a true phase transition while the other transition a. is not. 
The first part of this expectation may be justified by the observation that phase transitions 
occur in other models in which there are vastly more high-energy states (random graphs) than 
low-energy states (low-energy graphs) |31i] . However, like in other statistical models, the details 
may depend on the dimensionality of the graph in the low-energy regime. The second part of 
the expectation can be based on the fact that no such dramatic increase in the number of states 
is present among the regular and irregular random graphs. In addition to these uncertainties, 
it is also not known whether the quantum version of the model (one in which superpositions of 
classical graphs are allowed as quantum states) would exhibit the said phase transitions. 

IV. SIMULATIONS 

The main obstacle in evaluating the partition function of a graphity model is in enumerating 
the allowed graph configurations and in estimating their associated energies. Similar difficulties 
with other statistical systems have led to the development of various numerical approximation 
techniques. Notably among these are different types of Monte Carlo simulations [1^. Such 
techniques can also be used to address various questions relevant to a graphity model. For 
example, simulations can search for the lowest-energy graph for a given set of parameters N, 
vq, and r and couplings gv and qb- They can also provide estimates for the partition function 
and the expectation values of obervables, such as the energy, for these parameters and couphngs. 
This section describes results obtained through simulations of classical graphity models (29| . 

Given the expectation, justified in Sec. IIIH that states should be fo-regular at the lowest 
energies, all the simulations treat exactly fQ-regular graphs only. This greatly simplifies the 
simulations and eliminates gy from the list of parameters that can be explored. The coupling 
associated with cycles qb, thus becomes the only coupling in the basic Hamiltonian and can 
thus be normalized to unity, 

9B = 1. (46) 

Despite the restriction to fo-regular graphs, the scope of the simulations must be further 
restricted for computational reasons. In particular, since counting cycles requires an algorithm 
that runs in exponential time in the length of the cycles, simulations are practically restricted 
to small values of L^ax and thus of r. 

The results presented were obtained using the Metropolis algorithm. The moves used to 
alter the graph during the simulation are those shown in Figs. [T](a) and Fig. [UJ^b), selected at 
random at each time step of the simulation. For a wide range of temperatures, the Metropolis 
algorithm performs well and produces results quickly. However, the graphity Hamiltonian 
does in principle produce an energy landscape that has multiple local minima, and hence the 
Metropolis algorithm encounters the danger of getting trapped in one of those local minima 
at very low temperature and not exploring the full configuration space. In this work, this 
problem is avoided mainly by using relatively small N, avoiding extremely large value of (3, and 
repeating suspect simulations multiple times if needed. In some cases, the simulated tempering 
technique has also been tried, but all the results reported below rely only on the Metropolis 
method. 



15 




(a) (b) 

FIG. 2: Lowest-energy graphs in the basic model with = 36 and with (a) r = —2.5 and (b) r = 2.5. 
The results are obtained through simulation at /3 = 0.2. 

A. Basic Model 

When studying the basic model only, the general graphity partition function (fT8|) is restricted 
using 

z{G) = l, C{G)=(3E{G). (47) 

Within this framework, consider models defined by the parameters Vq = 3 and r = ±2.5. 
The first parameter vq restricts the simulation to cubic graphs. The low value of the second 
parameter r suggests that only short cycles contribute significantly to the energy of a graph. 
However, since the Hamiltonian depends on cycles of various lengths, the length cutoff L^nax 
must be set higher than the minimal cycle length. In practice, setting Lmax = 9 allows to 
compute a graph's energy at the percent level. 



1. Low- Energy Graphs 

One set of simulations can be performed with the aim of finding the lowest-energy graph 
configuration. To do this, the initial state of the simulation should be sufficiently different from 
the expected minimal energy graph; random graphs make for good initial states and can be 
generated efficiently [30|. The simulation should further be performed with N and f3 suitably 
chosen so that the system evolves rather quickly toward low-energy states, while still allowing 
some fluctuations to occur. Fluctuations ensure that the true, rather than a local, energy 
minimum is found. After some trial and error, the values = 36 and f3 = 0.2 are found to 
satisfy these criteria. The lowest-energy graphs found with these settings are shown in Fig. [21 

The two configurations shown are very distinct. When r is negative, the system effectively 
maximizes the number of 4-cycles in the graph. This leads to the graph in Fig. [2](a) containing 
six groups of six vertices connected in a chain. For positive r, the system allows some 3- and 
5-cycles, and this results in the treelike structure shown in Fig. [2]^b). Curiously, one of the 
leaves in Fig. EJ^b) is different from all the others. This is because for = 36 there can be no 
graph in which all the leaves have the most optimal five vertex structure. The position of the 
unique leaf on the tree is arbitrary, and thus the lowest-energy graph is degenerate. 
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FIG. 3: Average energy per node as function of in two series of simulations. Data points in the 
rising series are obtained through simulation with vq =, r = —2.5, and fixed /? = 0.17. Data points 
in the flat series correspond to /? scaling logarithmically with N. Error bars represent the standard 
error on the mean. The lines shown merely connect the data points and do not represent any fit. 



The general features of Fig. [2] persist for other values of A^, Vq, and r. The lowest-energy 
states of systems with negative r are always chains of some small components. The internal 
structure and the size of these components depend on the particular choices of parameters, 
but the overall one- dimensional chainlike structure does not. For positive r, the lowest-energy 
structure is always a tree. Again, the structure of the leaves varies with r. 



2. Scaling 

The scaling estimates of Sec. IIIIBI suggest that systems should become increasingly domi- 
nated by random graphs in the large A^ limit unless the inverse temperature j3 scales at least 
logarithmically with A^. This effect can be verified using Monte Carlo simulations. 

Consider performing two series of simulations of the r = —2.5 model. In the first series, 
simulations at different A^ are performed always using /5 = 0.17. In the second series, the 
temperatures are chosen according to (3 = jSolnN with /3o = 0.17/ In 30 so that the two series 
match for A^ = 30. In each simulation of each series, the expectation value (in the sense of the 
canonical ensemble) of the cycle energy per node, {E{G))/N, can be computed by extracting 
a number of representative configurations from the simulation and averaging them [28|. Since 
the lowest energy per node is independent of A^ for large N, this expectation value can be 
meaningfully compared across simulations with different A^. The results are shown in Fig. [21 

In both series, the average energy per node is negative for A^ = 30. This indicates that the 
system is dominated by configurations close to the minimal energy state (the lowest possible 
energy per node for r = —2.5 is eo = —12.2). In the series with fixed /?, however, {E)/N rises 
sharply with increasing A^. This reveals a transition to a regime of the parameter space where 
random graphs start dominating over low-energy graphs. The line does not reach {E)/N ~ 0, 
because for these still small values of A^ there are still a large number of graphs that have quite 
a large number of cycles but no longer resemble the chain structure. In other words, points 
with eo < {E)/N < represent graphs of the type omitted in the estimates of the partition 
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FIG. 4: Average energy per node for systems with r = —2.5 and equal to (left to right) 60, 120, 180, 
and 240. Error bars on the points (some too small to be seen) represent the standard error on the 
mean. The lines do not represent computed data and are only shown to aid the eye distinguish between 
the four series. 



function in Sec. IIIII 

The second series shown in Fig. |3] is the one where /3 scales with A^. In this case, the 
average energy per node remains very small even for large A^ indicating that the low-energy 
configurations continue to dominate. Note that the dip at A^ = 150, which seems to conflict 
with the slightly decreasing trend, is a statistical artifact - the point is actually within 2a of 
the value corresponding to A^ = 180. However, the dip does indicate that the simulations start 
becoming less reliable at high A^ and increasing f3 because of the ruggedness of the energy 
landscape. This unreliability is a systematic error that is not included in the error bars shown 
on the plot. 

The data in Fig. [3], as all numerical work, cannot prove the assertion that the required 
scaling of (3 must necessarily be logarithmic nor that the logarithmic scaling is sufficient for 
the average energy per node to remain close its minimum. This is particularly so because the 
explored range of A^ is small given that the expected dependence of the transition temperature 
on A^ is logarithmic. Nevertheless, the data show the Monte Carlo simulations of the full model 
to be consistent with the expected behavior derived in section UTTl and thus provide an important 
consistency check for those estimates. 



3. Transitions 

Apart from verifying the estimated scaling of j3, one can also study the low-energy to random 
graph transition in the standard manner by performing simulations at various (3 for fixed r and 
A^. The results are shown in Figs. [Hand [5] for systems with r = —2.5 and A^ = 60, 120, 180, 
and 240. 

The plot in Fig. H] shows how the average cycle energy per node decreases as (3 increases. 
The transition is quite smooth for A^ = 60 but becomes sharper for larger A^. This effect is 
mainly due to the fact that the difference in numbers of states in the low-energy and random 
regimes becomes more pronounced with large A^. Also, random graphs with A^ = 60 might 
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FIG. 5: Normalized specific heat for systems with r = —2.5 and equal to (left to right) 60, 120, 180, 
and 240. Error bars are computed using the bootstrap resampling method. The lines do not represent 
computed data. 



have quite a few short cycles thus lowering the average energy at low (3. The location of the 
transition temperature shifts slightly to larger j3 with increasing iV, which is consistent with 
the scaling of [3 found in Sec. IIIIBI 

Figure [5] shows the specific heat associated with the same simulations. The values plotted 
are actually specific heats divided by A^^ to make all four curves have roughly equal height. It 
should be noted that no special effort was made to find the true maxima of the specific heats. 
Doing that would require more computer time due to the "critical slowing down" phenomenon 



within the Monte Carlo simulation j28j. However, the location of these maxima can still be 



estimated from the plots: they occur roughly at /? = 0.13,0.16,0.17,0.18. These numbers can 
be compared with the estimates from Eq. fl43al) . Using eo = —12.2, one finds the expected 
transition inverse temperatures to be /? = 0.17,0.20,0.21, and 0.22. The estimates therefore 
consistently overshoot the real transition values by roughly 20%. While the agreement is not 
stunning, it does show that the estimate can be a useful guide for determining where the 
transition occurs. 



B. Matter Using Variance 

The general partition function of a graphity model, Eq. f[T7|) . includes a contribution z{G) 
from the matter degrees of freedom. While the form of this matter factor must depend on 
how the matter is implemented in the model, it is reasonable to assume that this factor has 
some general features. It may be reasonable to assume, for example, that z{G) is large for 
homogeneous graphs and small for graphs that contain many inhomogeneities. This observation 
suggests to study the effect of matter by choosing an ansatz for z{G). A simple one is 

\yiz{G) = -iiV{G), (48) 

with /i a real number and 

y{G) = Y: {El{a) - f\ (49) 
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(a) (b) (c) 

FIG. 6: Graphs with lowest criterion C for = 36 and r = +2.5. Results are obtained with /? = 0.2 
and the variance potential set to (a) fi = 0.8 (b) fj. = 3, and (c) fj, = 4. 

the variance among the cycle energies EB{a) assigned to each vertex a in the graph G. This 
ansatz imphes using the Boltzmann criterion C{G) as 

C{G) = (3E{G)+ fiV{G). (50) 

Since n now appears hke a potential in statistical ensembles, it can be called the variance 
potential. When it is positive, inhomogeneous graphs are suppressed in the partition function. 

The variance ansatz for z{G) is not intended to mimic a particular type of matter. Rather, 
its main advantage is that it can be implemented in a straightforward fashion in Monte Carlo 
simulations. Searches for the lowest-criterion graph can now be conducted similarly as in Sec. 
IIV A II for varying values of the variance potential fi. 

1. Cubic Graphs 

Consider augmenting the simulations from Sec. IIV A II (A^ = 36, Vq = 3, and r = ±2.5) by 
the matter ansatz with variance potential fi. The results for r = +2.5 are shown in Fig. O 
The graphs shown are more and more homogeneous for larger /x, and the ladder graph shown 
in Fig. M^c) is entirely vertex-transitive. Interestingly, this graph contains no cycles of length 
three in sharp contrast to the ground state graph obtained when fi = shown in Figure [2](b). 
Simulations of the system with r = —2.5 with large fi yield the same ladder graph. 

In terms of emergence of geometry, the results obtained show that an effective homogeneous 
one-dimensional space can be formed from the graphity model with r = —2.5. Moreover, for 
either sign of r but high enough /i, the emergent space can be the same. 

Since the simulations leading to Fig. E] are performed at low temperature, they can also give 
some indication as to what the next-to-lowest-energy graphs look like. Some of these graph 
excitations for the r = —2.5 simulation are shown in Fig. [3 These excitations quickly start 
looking like randomized graphs. Nonetheless, one can still recognize the background lowest- 
energy graph from their diagrams. 
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(a) (b) 
FIG. 7: Two graphs representing excitations of the r = —2.5 and high fi model. 




(a) 



(b) 



FIG. 8: Low-criterion quartic graphs, (a) is obtained with r 
produced with r = —2.5, P = 0.17, and /i = 3.8. 



-1.5, /3 = 0.6 and // = 32. (b) is 



2. Quartic Graphs 

Quartic graphs are graphs in which all vertices have degree four. Simulations with vq = 4 
can be performed in the same manner as for cubic case. One new property of quartic graphs 
(and generally of regular graphs of degree larger than 3) is that long cycles may visit a given 
vertex more than once. 

Quartic graphs are of interest because they can give rise to several types of latticelike con- 
figurations. In particular, the square lattice corresponds to flat 2D space and has four-sided 
plaquettes. Also, the diamond lattice corresponds to flat 3D space, and its shortest cycles have 
length six. 

Simulations searching for graphs with lowest criterion C {G) were performed for small values 
of r and large variance potential /i. Some of the results are shown in Fig. [Si Neither of 
the shown graphs is planar and this indicates that the emerging graphs might be interpreted 
as having dimension greater than one. Indeed both graphs in Fig. [8] may be thought of as 
two-dimensional tori. This is particularly convincing in the case of Fig. [8](b), where the two 
circumferences associated with a torus have lengths four and nine. 



3. Graphs with Larger Plaquettes 

Simulations with |r| much greater than 2.5 encounter at least three problems. First, as |r| 
increases, the evaluation of the graph energy becomes more time-consuming as longer cycles 
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must be taken into account. Second, in order to obtain results that can be generalized to the 
continuum limit, the number of vertices used in the simulation must be large enough so 
that in the high-/3 regime, the energy per node of the emerging graph is not skewed by the 
finite value of A^. And third, as |r| increases, the depth of the valleys in the energy landscape 
increases. All three of these issues make it more difficult for a simple desktop simulation to 
provide accurate results quickly. More sophisticated methods like simulated tempering should 
likely be used to make real progress. 

It is possible to perform some quick simulations with large |r|, for example r = —6, if L^ax 
is set very small, for example Lmax = 7. Such a low cutoff value for the cycle length means 
that the evaluation of the energy is not at all accurate. However, all three problems mentioned 
above become less serious. Simulations performed in this way with large f3 and large /i produce 
graphs with hexagonal plaquettes that resemble those graphs described in [l8[. While this is 
encouraging in the sense that it suggests that interesting lattices may appear in other regions 
of the graphity parameter space, it should be stressed again that such simulations with low 
Lmax cannot capture properties of the full graphity Hamiltonian. Proper simulations at large 
|r| must thus be considered beyond the scope of this initial numerical exploration. 

V. STATISTICAL MECHANICS OF A GRAPHITY MODEL WITH PARTICLES 

In the previous section, the effect of matter on graph configurations was modeled by adding 
a term proportional to the vertex energy variance to the Boltzmann criterion. Simulations 
showed that large values of the variance potential can indeed change the type of graph that 
minimizes the Boltzmann criterion ([T^ and thus change the structure of the expected graph 
at low temperatures. But while the variance prescription is straightforward to implement in 
a Monte Carlo simulation, it does not fully incorporate the effect of matter on a graph. It 
does not, for example, account for the fact that spaces with different geometries may allow for 
different densities of particle states. Thus, it would be desirable to supplement, or replace, the 
variance prescription with something more realistic. 

Since explicitly evaluating the partition function of a graphity model, with or without matter, 
is not feasible without approximations, one cannot, at least at this stage, learn very much by 
studying a particular form of matter coupling to graphs. This section therefore focusses on 
general properties of wavelike matter and discusses how they affect a dynamical geometry 
such as in a graphity model. The main principles leading to the argument that matter can 
significantly alter, and even determine, the shape of a dynamical geometry are outlined below 
in section |VA[ These principles are then applied to dynamical graphs in Sees. IV Bl and IV CI As 
in Sec. Illlt only the leading terms in A^ and (3 will be written in the estimates for the partition 
function components. 

A. Flexible Boxes 

Consider a set of A^ particles in a two-dimensional rectangular box with fixed area A in 
which the lengths and Ly of the two sides can vary as long as L^Ly = A. As the geometry 
of the box is not fixed, the partition function of the system includes a summation over the 
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configurations, denoted {Lx,Ly), tliat tlie box can take. It can be written as 



^ z(Lx,Ly)e 



(51) 



where Eb{Lx, Ly) is the energy associated with a given box, and z{Lx, Ly) describes the matter 
content. In the case that Eb{L^, Ly) = for all and Ly, the properties of the system, 
including its geometry, arc determined by the matter factors and the space of configurations C 
over which the sum is taken. 

The space of configurations is the set C of possible pairs of lengths {L^, Ly). For the present 
purpose, the precise form of this set will not be important, and only some of its general features 
will be assumed and used. First, the set C should contain certain extreme configurations. It 
should contain a configuration very close to {L.j.,Ly) = (a/A, a/A), i.e. a square box. It should 
also contain some configurations that are as long-and-narrow as possible. Assuming that the 
shortest allowed length is i, these extreme long-and-narrow configurations are {i,A/i) and 
{A/i,i). Second, the total number of configurations in the set C should scale as some power of 
A/i"^. This condition can be implemented if the lengths L^ and Ly are discretized. It should 
be noted that while the minimal length scale £ appears in several of the formulae below as 
a regularization parameter, the limit £ — > can be taken at the end without altering the 
conclusion of the main argument. 

For the matter factors, consider z{Lx, Ly) to be partition functions for a collection of N 
particles. Thus 



The summation in Zi is over the energy states available to a single particle in a box with side 
lengths Lx and Ly. These states may correspond to the standing waves that can form inside the 
two-dimensional box and, if the matter is quantum mechanical, their energies may be associated 
with inverse wavelengths. For nonrelativistic particles with mass m, the energy states El^^l^ 
could be given by 



with Hx and Uy the usual integer quantum numbers. 

It is instructive to look at the lowest-energy particle states available in the extreme box 
configurations. In the large area limit these are 




(52) 



with 




(53) 




(54) 




(55) 



The discrepancy between these two energies may be very large when i is small. Minimal particle 
energies in nonextreme configurations lie between the two shown values. 
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With this particular choice of particle energy spectrum, the matter partition functions be- 
come 




While exact expressions for (l56l) are not known, they can be approximated by the number of 
states for which f3 times the particle energy is less than unity^. Such a prescription should at 
least provide the leading behavior of these factors. One finds 



In the first line, it was assumed that /3/(8m£^) > 1. 

While the leading behavior with A is the same in both cases, the contribution from the 
square box can be larger than those from the long-and-thin boxes when i and m are small and 
(3 is large. In such cases, the dominance of the square configuration will be further amplified 
when Zi{Lx, Ly) is raised to the power of in fl52p . Of particular interest is taking N to scale 
with some power of A, i.e. ~ cA""; the special case a = 1 corresponds to considering a set of 
boxes of various areas A with fixed particle density. In the case with ~ A, the matter factor 
z{Lx, Ly) will be exponentially larger for the square configuration than for any of the other 
long-and-thin configurations. Moreover, since the number of nonsquare configuration scales 
as a power of A/i^, the contribution of the square configuration will dominate the partition 
function for this toy model in the limit A —>■ oo. Similar scaling arguments can be used to show 
that expectation values of observables and their fluctuations will be determined primarily by 
their values on the square configuration. 

This example is an example of how matter degrees of freedom can determine the shape of 
a dynamical geometry. The result is quite striking - out of the large set of configurations C, 
matter can select a special configuration and effectively force all observable physics to be based 
on this special geometry. 



The discussion of particles in a flexible box may be adapted to the case where the trap- 
ping container is a graph. The purpose of this section is to estimate the matter factors z{G) 
associated with different classes of graphs in the same spirit as zi{Lx,Ly) and z{Lx,Ly) were 
estimated above. In all cases, the number of particles, denoted N above, will be assumed to 



In the special case that — Ly, zi may perhaps be better estimated using an integral expression 




(57) 



B. Configurations 




n dn 



with f3 being a rescaled version of /3. This approach, however, does not significantly differ with the other 
scheme. 
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be the same as the number of vertices present in the graph. Thus one of the key assumptions 
in the following analysis is that the density of particles is constant across graphs with different 
number of vertices A^. 



1. Homogeneous Graphs 

Homogeneous graphs are the easiest ones to deal with since their large-scale properties 
resemble lattices. Because the graphity models all have a finite number of vertices, these 
homogeneous lattices are closed geometries and, as discussed in Sec. IIVB| represent two- 
dimensional tori and similar compact structures in other dimensions. 

Consider first two-dimensional tori. They can be described by two lengths Ci and C2 rep- 
resenting the two circumferences. The problem of evaluating z{G) for such graphs is in direct 
correspondence to evaluating z{Lj;, Ly) using fl56l) and fl52l) above. One modification that needs 
to be made lies in replacing the area by the number of vertices, 

A = Nf- (58) 

the factor £^ is required for dimensional reasons. Also, since each particle's wavefunction as 
well as its derivative should be continuous on the whole torus, the summation over particle 
energy states should be restricted to even integers. This restriction only rescales the estimates 
( !57|) by some factor of order one and thus the estimates made in the present case are equivalent 

to dnzD. 

To obtain the full contribution of the lattice graphs to the partition function, the estimates 
for z{G) must be combined with the nonmatter components of (l38l) . For a torus with two 
circumferences almost equal, Ci ~ C2 ~ iVN, one obtains 

In ^l;'^^,^ ^ N\nN-(3N{gv + eH) + N\n (59) 



For a long-and-thin torus, Ci <^ C2 or C2 ^ Ci, one has 



In ^l:'c.«c. -N\nN-(3N(^gy + aeu + ^) + \n In ■ 



(60) 



In both cases, note that the factor = {Ni'^)'^ appearing after raising zi to the power of 
is cancelled by the A^! in the denominator of (1521) . In the equation for the long-and-thin torus, 
the factor a appearing in front of e/^ is a real number of order one that reflects the fact that 
the energy per node associated with the two types of tori may be different. In particular, since 
the long-and-thin torus contains winding cycles that are not present in the isotropic torus [3] , 
this factor should be expected to be greater than one. 

Next, consider homogeneous graphs that represent generalizations of tori in higher dimen- 
sions d. For simplicity, consider the isotropic configurations in which all the circumferences 
C ~ iN^^'^. The generalization of fl56|) for this case can be written as 



(61) 
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with n being a (i-dimensional vector of quantum numbers. Using the same approximation as 
before, that the sum is approximately equal to the number of states such that the exponent is 
less than unity, and assuming that P{8mi'^)~^N~'^^'^ < 1, one obtains^ 

z^C ~ m'/') ~ N l^^j , (62) 

i.e. the same dependence on as given by the expression (1571) . This is quite important: 
despite the number of states on a (i-dimensional space up to some level \n\ growing as \n\'^, this 
growth is largely balanced by the smaller lengths L ~ iN^^'^, resulting in an overall Zi that 
does not depend very strongly on d. After using this to evaluate the full contribution of the 
c/-dimensional torus to Zh, one therefore finds an expression that has almost the same form as 
(l59l) . namely 

1 8m£^ 

lnZ^ ,,^^i/rf~ NlnN - PN {gv + a^en) + -Ndln^. (63) 

The most important difference between this expression and (159|) lies in the different energies, 
encoded in the real numbers a^, associated by the graphity Hamiltonian to the different tori. 

Estimating the contribution of nonisotropic (i-dimensional tori to Zh can be done in a similar 
fashion. As in the case of two-dimensional tori, such nonisotropic structures would contribute 
terms which, like ( l60l) . would contain factors exp {—j3N/{8mi'^)). All the different terms will 
not be written out here as they will turn out to not to play an important role. 



2. Low-Energy Graphs 

Simulations show that the lowest-energy graphs are one- dimensional chains made up of 
subgraphs with a small number of vertices. Since such graphs are not homogeneous lattices, 
one has to make further assumptions in order to guess what their matter factors should be. 
Suppose that these chain graphs can support particle states as standing waves both along the 
one-dimensional chains and also within the individual subgraphs. Then, since the substructures 
only have a finite number of vertices, standing waves supported within the subgraphs will have 
wavelengths no longer than a few i. Thus, these standing waves will have a large energy; the 
situation is not dissimilar than for the long-and-thin boxes or tori. 

This reasoning suggests to write the contribution Zl with a factor that penalizes chain 
graphs for high-energy particle excitations, 

InZLr^ NlnN -PN (^gv + eo + ^^y (64) 

Here, ttil is an effective mass corresponding to the particles on the low-energy graph. It is 
thought to absorb all the factors that are missed in this simple estimation. What will be 
important is that this effective mass should not scale significantly with A^ and should always 
be positive. 



^ If the sum were approximated using the integral approximation in footnote 1, this would have an additional 
factor depending on d. This discrepancy will not be important for the argument below. 
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3. Random Graphs 

As in the case of low-energy graphs, it is unclear how to estimate the matter factor z{G) 
associated with random graphs. In general, the term for fo-random graphs can be written as 

In Zr,,, ~ ]^Nvo In iV - jSNgy + In z{Gr) . (65) 



Since the average distance between two vertices in a random graph is small [27[, the largest 
wavelengths associated with particle states should be on the order of ^. This would suggest to 
write this factor similarly as in flMl) . for example, 

with some effective positive mass tur. Such a matter factor would be effectively penalizing 
random graphs for their tight connectivity. 

A discussion of random nonregular graphs would follow the same lines. It will not be 
necessary, however, to include it here. 



C. Scaling 

Armed with estimates for the important components of Z, one can repeat the analysis of 
Sec. IIIIBI for the graphity model with matter. However, given that it was found previously 
that basic intuition can be extracted quite easily just by comparing the leading order behavior 
of the relevant terms, only this shortcut will be followed here. 

Because of the additional factors due to matter, the low-energy behavior of the system can 
be quite different than before. Compare, for instance, the partition functions (!59|) and (!60|) 
associated with homogeneous graphs. The contribution due to the isotropic torus can be much 
larger than the contribution from long-and-thin torus if (recall that aen < < 0) 

1 1 , 8m£2 

^-<"^- + 8;;^ + 2^^"-r- 

The last term on the right hand side becomes negligible if £^ is small and (3 large (as it is 
expected to be if it scales logarithmically with A^). The remaining inequality may be satisfied 
by judiciously choosing the parameters m, and the coupling qr appearing inside the energy 
en- Thus, matter can select the isotropic graph in the low-temperature regime even though the 
long-and-thin configuration contains more cycles. Next, comparing quantities (1631) for various 
(i, one can find the dimension d of the most optimal torus. This dimension d will be primarily 
determined by the number of short cycles, i.e. by the numbers ad- 

Suppose that the optimal dimension is d = 2 so that the dominant homogeneous term is 
given by fl59l) . Compare this to the contribution flM|) from the low-energy chainlike graphs. 
The homogeneous graph can be dominant if the couplings are chosen such that (recall that in 
general eo < e^/ < 0) 

1 1 , 8m£2 ^ ^ 
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If this is satisfied, then low-energy behavior of the graph system with matter is drastically 
different from what was found in Sec. IIIII Also, as the matter contribution to the partition 
function becomes important, the cycle counting becomes less so. This implies that the un- 
derlying graph structure, while preserving the large-scale properties of the emergent geometry, 
can perhaps become more fluid at small scales than the figures of homogeneous lattices might 
imply. 

In the high-energy regime, just as before, random graphs can dominate above some transition 
temperature. The transition temperature to this regime should again be proportional to In N. 
The estimate of the coefficient in front of this factor should depend on the effective masses 
m and mn of the particles in the high and low-temperature regimes. Since the matter factor 
associated with random graphs are not known in detail, this transition is not discussed further 
here. 



VI. DISCUSSION 

Summarizing, this paper focused on the classical statistical mechanics of graphity models 
and produced two types of results. First, results from Monte Carlo simulations show that 
neither r > or r < basic models produce ground state graphs that are latticelike. The more 
interesting case with negative r can at best produce an extended graph that is made up of 
smaller units chained together in a circle. Thus it appears that the graphity Hamiltonian that 
depends on vertex degrees and counts closed cycles must be supplemented with some additional 
conditions. Indeed, ground states searched for when the original graphity model was introduced 



were implicitly assumed to be homogeneous or almost homogeneous llSH . Simulations of 
extended models implementing such homogeneity conditions in turn reveal that the intuitions 
in the early works [l7, 3] are correct and that emergent graphs in the high inverse temperature 



regime can really describe extended geometries. 

Second, the partition functions for the models with and without matter are estimated ana- 
lytically to better understand the numerical results. In the context of the basic model without 
matter, such estimates, while very rough, can nonetheless provide useful information about the 
behavior of the graphity models in the continuum N ^ oo limit. If this limit is defined through 
a family of models with constant couplings, the transition inverse temperature f3 can be shown 
to scale logarithmically with A^. Therefore, the transition occurs at zero temperature for large 
A^. (The transition temperature could be finite if, for example, the couplings were taken to run 
logarithmically with A^ instead.) In the context of the model including matter, estimates of 
the extended partition function explain how matter degrees of freedom can alter the low and 
high-energy behavior of a graphity model. In particular, they show that matter can enhance 
the importance of homogeneous graph configurations in the partition function to the extent 
that these graphs become dominant. Unfortunately, it was not possible to test the particular 
form of matter used in these estimates in numerical simulations. 

The problem with nonextended configurations in the large inverse temperature regime in the 
graphity models is reminiscent of similar findings in the context of Euclidean dynamical trian- 
gulations j^. In that context, the problem was noticed to become milder when the triangulation 



was coupled to gauge fields [21|, |22|, but the main solution has been based on restricting the 



configuration space to those triangulations respecting certain criteria motivated by causality 
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jsl and the presence of global spatial hypersurfaces 3^. In graphity models, such a restriction 
of the configurations space cannot be naturally introduced. Instead, the proposal to bypass the 
problem of nonextended geometries is to make the configuration space larger to include matter 
degrees of freedom. This is close in spirit to the mentioned work on dynamical triangulations 
with matter fields 21, 2^', the very interesting question of whether the mechanism described 



here might also regularize Euclidean dynamical triangulations deserves some further thought. 

In the end, what can be learned from simple models based on graphs? At the very least, 
they should be seen as concrete examples of setups that do not require a continuous mani- 
fold as an ingredient but rather are able to produce one (or an approximation to one) in an 
appropriate regime. In particular, they can produce flat emergent geometries in the form of 
tori. Given this, one can start asking if some graphity model can actually reproduce the known 
gravitational dynamics in three plus one dimensions in some region of its parameter space, or 
if the discreteness may have observable consequences for particle physics or cosmology (works 



studying related issues in similar models include [33|, l3J].) These are issues that are worth 
exploring further and now bring graphity models, which begin at a slightly speculative starting 
point, possibly close to being compared with experiment. 
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